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Abstract 



' The spectral decomposition is given for the N-qubit BeU operators with two observables per qubit. 

It is found that the eigenstates (when non-degenerate) are N-qubit GHZ states even for those operators 
that do not allow the maximal violation of the corresponding inequality. We present two applications 
QQ ' of this analysis. In particular, we discuss the existence of pure entangled states that do not violate any 



Mermin-Klyshko inequality for A'' > 3. 



1 Introduction 



^ • There is a huge literature about the topic of Bell's inequalities. For the sake of this introduction, we briefly 
' recall some well-known ideas, using the Clauser-Horne-Shimony-Holt (CHSH) inequality Q. For each choice 
, of four numbers ai, a'j^, a2, a'j € { — 1,+!}, the quantity S — h{o-2 + a2)'^i + 5('^2 ~ fl2)'^i can take only 
the values -1-1 or -1. Therefore, if the four numbers are considered as realization of random variables, the 
expectation of S will certainly depend on the distribution of the variables but must satisfy |£'(>5')| < 1. 
' These are "trivial mathematics". But if one turns to quantum mechanics (QM), then this inequality can be 
5t 1 violated. Indeed, consider that a is +1 if, as a result of a measurement, a spin is fomid along the direction 
-|-a, and a is —1 if the spin is found along the direction —a. This is achieved by replacing a by the operator 
a • (T = (Tq. Using this prescription for ai, a[, a2 and flj we find that S is the expectation value of the "Bell 
operator" [|| 

B2(ai,a'i,a2,a^) = i (ctq^ + Cr„^) ffai + ^ (era, - (Ta'J «) CTq; (1) 

where ai , a'j^ , a2 , a2 are unit vectors that will be referred to as the "parameters" of the Bell operator. It is 
well-known that for some choices of the parameters the highest eigenvalue of B2 can be higher than 1: there 
exist some states j^*) that violate the inequality |(-B2)*| = < 1- "Trivial mathematics" fail whenever 

[(Tdfe , CTaj^j 7^ 0, that is whenever at 7^ i^k- ^^^^ case the "random variables" ak and a'j, cannot have 
simultaneously a precise value. 
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The first Bell's inequalities were derived for two two- level systems (hereafter referred to as "qubits"). Gener- 
alizations of the Bell's inequalities have been proposed along the following lines: (i) bipartite inequalities for 
two n-level quantum systems (ii) bipartite inequalities using more than two parameters per system ai, 
a'l, a'(,... Q; (iii) multipartite inequalities, that is, inequalities involving more than two quantum systems 

In this paper, we consider inequalities involving an arbitrary number of qubits using two observables per 
qubit (the observables are obviously dichotomic). This family of inequalities has been studied in great detail 
independently by Werner and Wolf jsj and by Zukowski and Brukner |^ . Our present contribution consists 
in exhibiting explicitly the spectral decomposition of the Bell's operators (section 2). Section 3 presents 
applications of this result. 



2 Spectral decomposition of Bell operators 
2.1 Inequalities for two observables 

Consider a quantum system composed of n qubits, that is a system described by the Hilbert space (C^)*^". 
For each qubit fc, we define two observables Ak(0) — aa^. and Afc(l) — aa'^, with at and aj^ two vectors on 
the unit sphere. The set {ai, a'j^, ...a„, a^} of the 2n unit vectors is written a. Up to normalization, any 
n-qubit Bell inequality]^ can be written as (;B„(a)) < 1 for a given Bell operator Bn- The form of the Bell 
operator is a polynomial 



n 



Bnia) - Pi^)<S)Msk)- (2) 

se{oa}" k=i 

The coefficients /3(s) are rather arbitrary, provided that < 1 is satisfied for all product states. 
Of course, not every polynomial of the form (H) defines a good inequality; in the worst cases, e.g. when the 
polynomial is simply (Tai ^ ■■■ ^ cra„, one will find < 1 for all states. The complete classification of all 
inequalities is due to Werner and Wolf Q . A special role is played by the Mermin-Klyshko (MK) inequalities 
^ , whose corresponding Bell operator is defined recursively as 

B„(a) = Bn ^ ^{cra„ + <Ja'J «) Bn-1 + ^{da^ - <Ja'J «) B'^.^ (3) 

where is obtained from S„ by exchanging all the a^ and a'j,. In particular, B2 is given by the CHSH 
inequality; B3 is the operator that corresponds to the so-called Mermin's inequality Q. For MK inequalities, 
the violation allowed by QM is (Bn) = 2*^"^^^/^; no other inequality with two observables per qubit can reach 
such a violation [^. More results on these operators are given in Appendices A and B. 

2.2 Spectral decomposition: statement of the theorem 

We want to characterize the eigenvectors and eigenvalues of the Bell operator Bn defined in (^) , for a given set 
of 2n unit vectors a. For this purpose, we can suppose without any loss of generality that all the unit vectors 
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in a lie in the (x, y) plane; physically, this amounts to say that the axes x, y, z can be defined independently 
for each qubit. We will show the following 

Theorem 1 Let Bn given by (Qj, with at = cosa^^ex + cosafeBy and aj^. = cosQ!'j.ex + cosaj^ey for all 
k = Let |0), resp. |1), be the eigenvector of for the eigenvalue +1, resp. —1. Finally, let 

n = {uji, ...,U!n) G {0, 1}" be a configuration of n zeros or ones, and — {lji, w„) with (Iifc = 1 — cjfe the 
complementary configuration. Then: 

L The 2" n-qubit GHZ states, labeled by the configurations fl, defined by: 

l^niOn)) = ^ (e^'-l^^) + 1^)) (4) 
form a basis of eigenvectors of Bn for some 9q ~ (^ni^)- 
2. The parameter 9fi — OQ{a.) and the eigenvalue Afj = Ao(a) are calculated from a complex number fn{a.): 
«/ arg/n e [0,7r[ then: flf^ = - arg /o - tt , \n = -\fn\; 

(5) 

if arg/n e [tt, 2tt[ then: 9n ^ - arg/o , • 
The complex number /f2(a) is obtained as follows: take Bn, and for all k = I, ...,n replace the operator 
Ua^. by the complex number e*"'' if ujk = in fJ, or by the complex number 6^*"*= if uJk — 1 in Q; and 
the analog replacement for Ua'^ ■ 

About statement 1: It was noticed in (V, D) that GHZ states are the states that maximally violate 
any inequality. This is an immediate corollary of statement 1, since for any matrix M it holds that 
maxi, {v\M\v) / {v\v) is the maximal eigenvalue, obtained if and only if v is the eigenvector associated to 
that eigenvalue. 

About statement 2: The definition of fni^.) may seem cumbersome, but a single example will clarify it. Take 
the CHSH operator is given by @). Then /oo = ^e'"^ + e^°''^) e'"^ + ^(e*"^ - e"^)e"i; /oi is obtained by 
replacing e*"^ and e'"2 by their conjugates; /lo = /gx and /n — /qq, with z* the conjugate of z. 
The proof of the theorem is given in two steps. In the first step, we take advantage of a remarkable symmetry 
of the Bell operators to guess the basis (||); in the second step, direct calculation allows to get the announced 
exphcit results (||). 

2.3 First step 

For a given rotation matrix R G SO (3), it is well-known that one can find U E SU (2) such that U a.-a U^^ = 
(i?a) • a. In particular, one can find U such that UuaU^^ = u-a — ^Oa- Since we are considering that 
all the parameters a lie in the plane [x,y), the rotation that brings a on —a is a rotation by tt around 
the z-axis, so that the corresponding unitary operation is C/ ~ CTz (equality up to an arbitrary phase). We 
introduce the notation L/[fc] = 1 ® ... (g) cr^ ... ® 1, where the rotation is applied on the fc-th qubit. Note 
that [?7[/c] , ?7[/]] = for all fc and I. Since Bn is a sum of terms like aa^aa',^...aa„, we have manifestly: 

U{k\BnU{k\-^ = -Bn Vfce{l,...,n} (6) 
U\k\V^\BnV\k\-^\J\lY^ ^ Bn Vfc,/e {!,... (7) 
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These conditions depend critically on the assumption that the Bell operator is dichotomic: if we had three 
or more vectors for a qubit, we could not ensure that they lie in a plane. Condition says that Bn and 
—Bn are linked by a unitary operation, whence the following: 

Lemma 1 If X is an eigenvalue of Bn associated to lip), then —A is also an eigenvalue of Bn- The vector 
U[k]\ip) is eigenvector of Bn for the eigenvalue —A, for all k. 

The symmetries and (|^) of Bn suggest to look for vectors satisfying 

UM^)) ±1*) Vfc < n; (8) 
U[k]U[l]\^) = e'-'^'l*) yk,l<n (9) 

as good candidates for the eigenstates of Bn — they would be the unique candidates if none of the eigenvalues 
of i3„ were degenerate, but this is generally not the case (see Appendix B). Let |0) (resp. |1)) be the eigenstate 
of (Tz for the eigenvalue +1 (resp. -1). We decompose the n-qubit state j'l') on the basis of the product states 
of |0)s and |l)s: l^*) = Ene{o i}" We use condition (^ first: 

U[k]U[im = ^(-ir-^+^'cnir!) e'^'"Y,cn\n)yk,l ^ 
n o 

^ [e'^'^'^+'^i) ^e''""]cn = yne{0,ir and yk,l <n. (10) 
Now suppose cji =/= 0: this implies, modulo 2tt: 

{0 if ujk = 
(11) 
TT if UJk 7^ 

That is, the choice of for which cn 7^ determines completely the sequence of the jm- Now, it is 
evident from (|ll|) that only fl gives exactly the same sequence as fl. Thus ( |lO[ ) means that once we have 
chosen ft for which cq ^ 0, then cji' = for all il' ^ fl, 0,. We turn now to condition (||), that, with 
t^[fc]|*>-Eo(-ir'co|f^),reads: 

i^lUm = ^(-l)"^|cnp OVA:. (12) 
n 

But we have proved just above that the states we are interested in are such that only cq and Cq can be different 
from zero. Thus (|l^ ) becomes (—1)'^'= (jcop — IcjiP) = 0, that is \cn\ — |cq|. We have then proved that a 
n-qubit state satisfies both (||) and (||) if and only if it is of the form (H) for a given ft G {0, 1}". Thus we have 
2" states, each labeled by one configuration fl. The orthogonality requirement {'i'n{Sn)\'^n'{0n')) = ^n.n' is 
trivial but for fl' = 17: in this case, we must require 6^2 = tt — This concludes the first step of the proof. 
Just two remarks before turning to the second step: 

Remark 1: The state built on f2 is entirely determined by the state built on f2 through 

\^n{On)) ^ 1*0(^0 + ^)) ^ U[k]\^'n{en)) . (13) 

The first equality follows from the requirement 0^ — n ~ 6^2 by extracting 9^ as a global phase. As for the 
second equality: 

U[k]\^n{On)) = ^ {e^'-i-ir-in) + (-1)^1^)) = 
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= ^J' il^) - e-^'-m ^ j= [e^(--'-^\n) + \n)) . 

Thus, another recipe to build a basis of states of the form (^) is the following: (i) for all such as (say) 
wi — 0, choose 9q and build the state |5'(ri, 6*0)) ; (ii) apply U[k] to each of these states (or change 9fi to 
dn +71') to complete the set. 

Remark 2: The two states built from Hq = (0, ...,0), namely (e*^|0...0) ± |1. ..!)), are the only ones for 
which the "fki are independent of k and I — and are actually 0, with our choice of phases U[k] — az[k]. 

2.4 Second step 

We must show that the states of the form (jj) form a basis of eigenstates of Bn (a) , that is 

Bnl^fniOn)) = Xn\^n{On)) (14) 

with 9n — On{a) and Ao = Ao(a). Actually, we have to solve ( |l4| ) only for the state built on ilo = (0, 0) 

l^noiOn.)) ^ \^i0)) = ^ (e^'|0...0) + . 

This is so, because one can exchange |0) and |1) by applying a unitary operation, here a^. Therefore by 
application of to the suitable qubits we can always transform any l^'si) into l^'sio). Once we have the 
results for f^o, the results for Q follow by taking the qubits k for which cj/j = 1 in 17, and replacing aa^. by 
CajCafcCTx, that is, replacing {ak)y = sinafe by — sinafe, and the same for Caj^. So the eigenvalue problem is 
reduced to the problem of finding A and 6 satisfying 

S„(e'«|0...0) + 11...1)) = A (e^«|0...0) + 11...1)). (15) 

Consider now one of the terms in (||), say (8" ... (Ja„- a standard calculation gives 

'I0...0) + I1...1)) = e'^([|e'"'=)|l...l) + (]^e-'"'=)|0...0). 



CTq, (X) ... » (7a„ e 



k 



Consequently the eigenvalue problem (|T^ gives 

B„|f(0)) = i= (e^''/|l...l) + r|0...0)) = M'fie)) ^ e'Of = A (16) 

with /(a) = ( X^s /'(*) rifc ^"'°^^'°^) — of notation, we wrote afc(O) for ak and afc(l) for a'j,. The 

proof of Theorem ^ is virtually concluded. The solution (^ follows by settling a matter of convention, since 
condition ( p^ ) can be written as |/|e*(^+^''s(^)) = A or as |/|e*'^+''''s(/)+^) — —A; in other words, a convention 
on 9 fixes the sign of A — this is nothing but the manifestation of (|l^). We choose as a convention that 
0o G [0, n[ for all Q; this convention is consistent with 6^ = tt — ^o- 



3 Applications and perspectives 
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3.1 On some non-maximally entangled states 

In this section we study the violation of MK inequahties for a family of N-qubit states that clearly exhibit 
N-qubit entanglement. These states are 

\i1^n{4>)) = cos0|O^) + sin0|l^), (17) 

where we adopt the notation |0^) = |0...0); by convention, we choose cos0 > sin0 > i.e. (j) £ [0, j]. 
In the case N = 2, using Schmidt's decomposition every pure state can be written in the form = 
cos(?!>|OO) + sin0|ll). It is well-known that the CHSH inequality is violated by all pure entangled states in 
fact, using Horodeckis' theorem one can calculate explicitly that S2 — maxQ(_B2(a.))^,(0) = \/l + sin^ 20, 
which is bigger than 1 unless (j) — 0. It is interesting to re-derive this result starting from the spectral 
decomposition of B2. Consider the Beh states: |<I>±) = ^TjdOO) ± |11)), = ^TjdOl) ± |10)). |^(0)) is a 
linear combination of I^J) and 1$^). If we take the unit vectors a in the {x,y) plane, then |$+) and \^^) 
must be associated to opposite eigenvalues. But |(f>+) ~ |$+) and |$+) — I'i'^): therefore, by taking the 
unitary vectors in the (y, z) plane, we can construct B2 as 

B2 = Ai ~ P^-) + X2 (P^+ - P^-) = Ai (P^+ -P^+) + X2 (P^- -P^-)- 

This way, the two vectors that have a non-zero overlap with |V-'(0)) are associated to the positive eigenvalues. 
The calculation of S2 is not difficult, using the fact that Af + A2 = 2 (see Lemma ^ of Appendix A) and the 
standard maximization 

max(^ cos x + B sin X) = VA^ + . (18) 

X 

We find indeed Horodecki's value. Thus, to obtain this maximal violation of CHSH we took advantage 
of the possibility of choosing the Bell states that are orthogonal to \ip{(f>)) as the states associated to the 
negative eigenvalues of B2. Now, this is precisely a characteristic of two-qubit maximally entangled states 
that does not generalize to three or more qubits. In fact, it is well-known and easily verified that N-qubit 
GHZ states take the form -^(|0^) + |1^)) only in one basis (up to trivial relabelling). Therefore, for N > 2, 
if we build Bn such that |GHZ+) = "^(|0^) + |1^)) is associated to the eigenvalue A, then necessarily 
|GHZ_) = ^(|0") - |1^)) will be associated to -A. 

We now consider Sn — maxa(i?Ar(a))^j^(0), with Bn a MK Bell operator, for arbitrary N > 3. Using 
(l^|Bjv|l^) = (-1)^(0^ \Bn\0^) we have 

(P^(a))^0„(^) = /Ar((/))(0^|Sjv(a)|0^) + sin20Re((l^|i?Ar(a)|O^)) (19) 

where fN{4>) — [cos^ (j) + (— l)^sin^0], that is 1 for N even and cos2(/) for N odd. The maximization of 
( |l9|) over all possible choices of a is not evident for the following reason. We know that there are sets a that 
saturate the bound Rc((l^|PAr|0^)) = 2^; but for these we find (0^|PAr|0^) = 0. Similarly, the sets 
a that saturate the bound (0^|i37v|0^) = 1 give Re((l^|i3Ar|0^)) = 0. Let us try to guess the maximum 
of ( [1^ ) using the insight provided by the spectral decomposition of B]y discussed in section 2 above. A 
natural first guess would be Bjv = 2~^(Pqjjz ~ ^GHZ )' ^^^^ R.e((l^|i?Ar|0^)) = 2~2~. This choice 
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Figure 1: Ranges of (f) for which \iI)m{4>)) violates (brighter) and does not violate (darker) the MK inequality. 
For N > 5, the result is conjectured, see (|2C 



gives 5*1^ (0) = 22 sin 2(j). Numerical evidence suggests that this is indeed the maximum of (|T^) whenever 
'S'Ar(</') > 1- This provides a criterion for violation of MK inequalities: 

, . iV = 3,4,5: numerically verified . 
Spf{(l)) > 1 ^ sin20 > 2 — — I I . (20) 

iV > 5 : conjectured 

Thus there exist pure entangled states that do not violate the MK inequality. Let's define 0jv as the value of 

TV — 1 

(f> at which |-0Ar((/))) ceases to violate the MK inequality: we have (f>2 — 0, and sin2(j)]y = 2 ~ for > 3, 
within the validity of (pO|). There is a clear discontinuity in the behaviour of ^at between N = 2 and iV = 3, 
as illustrated in figure |l|: (pN jumps from for = 2 to for iV = 3, then starts decreasing again for higher 
A^. This analysis suggests that the MK inequalities, and more generally the family of Bell's inequalities with 
two observables per qubit, may not be the "natural" generalization of the CHSH inequality to more than 
two qubits. Whether a more suitable inequality exists is an open question. 

To conclude this section, let us briefly come back to the problem of maximizing (|l9|). The guess S^j^\(l)) = 
2T~sin20 cannot be correct for all 0: in fact, in the limit of small i/i's, |'0Ar(0)) approaches |0^) and 
therefore 5'(0) should converge to 1. To avoid this problem, we modify slightly our guess to have 

S'j^\(j)) - max [2^ sin 20, /at (0)] . (21) 

Obviously Sn > since wc know that S*^"* can be reached. Numerical estimates for N = 3,4,5 prove 

that this guess is extremely good. Actually, for A^ = 4, 5 we found SNi(t>) = s'"^\4') for all if) to within the 
accuracy of the calculation; for A^ = 3, the same holds for all (p but a small range of values around cj) defined 
by 2 sin 20 = cos 20 (the maximal difference is S3 — S^^'' ~ 0.002). 
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3.2 Bounds on the violation of Mermin's inequality 

As a second application, we show how the knowledge of the spectral decomposition of Bell operators provides 
bounds to estimate the violation of Mermin's inequality for any three-qubit state p. Some results are similar 
to those found independently by Zukowski and Brukner j^. 

We consider a three-qubit Bell operator 63 (a) . According to Theorem |l| one can always find a basis such that 
its eight eigenstates are |*i,8) = l^-ooolO, tt)) = ^(|000) ± 1*2,7) = |*ooi(0,^)) = ^(|001) ± |110)), 

1*3.6) = |^'oio(0,7r)) = ;^(|010) ± |101)), 1*4,5) = 1*011(0, tt)) = ;^(|011) ± |100)). Note that the four 
angles 9^ that are unconstrained can be chosen to be without loss of generality, since this amounts to a 
redefinition of the global phases of |0)^, |0)^ etc. Consequently in this basis 

B3 = Ai (Pi - Pg) + A2 (P2 - Pi) + A3 (P3 - Pe) + A4 (P4 - P5) = 

^ ^^++++'^xxx + M-+H — '^xyy + M — I h^yxy + M h+'^yyx 

with a^xx = <Jx ® cFx ® cTx etc., and with ^siS2S3Si — \{si\i + S2A2 + S3A3 + S4A4). For a given three-qubit 
state p 

Tr(S3 p) = p+++ + txxx + M — hH — ^xyy + P- — \ ^tyxy + P ^ + tyyx (22) 

with the standard notation t^xx — Tr(/3 <Txxx) etc. Our final purpose is to estimate Sp = max^ Tr(yB3 p) for any 
p. If p is given, one must find both the good eigenvectors and the good eigenvalues of B3. The optimization 
of the eigenvalues is performed by varying the parameters p; we discuss it in the next paragraph for the 
Mermin's operator P3. To optimize the eigenvectors means to define the axes x and y for each qubit. Note 
that when the basis of eigenvectors is optimized only four number tijk will come into play, thus sharpening 
the condition obtained by Zukowski and Brukner j^] that involved eight of these numbers. 
While the system of eigenvectors is the same for all Bell operators of the form (||) , the eigenvalues and their 
properties obviously depend on the operator that is considered. We restrict our discussion to the Mermin 
operator P3 given by (^). It can then be shown that the eigenvalues must satisfy 

Tr(P|) = 8 (23) 

(see Lemma || in Appendix A). This leads to p'^^^^ + + p'^^ ^ -I- p^ |__|_ = 1. Therefore, we can 

let = cos a cos /3, p = cosasin/3, p | |_ — sin a cos 7 and p ^+ = sin a cos 7, and maximize 

over a, /3 and 7. By using thrice the maximization {uM we find 



Sp < Sp — niax yt'^xx ~^ ^xyy + ^yxy "I" ^yyx ■ (24) 

This bound would in fact be exact if there were no constraint on the eigenvalues other than (^). However, 
starting from the eigenvalues as they are given in statement 2 of Theorem |^, one finds by inspection that 
the eigenvalues are bound to fulfill some other conditions, like ('^3+-^^i^^^)(^3^+-''4 ^) _ 2 sin^(Q:i — a[). To see 
that such a condition is indeed an additional constraint, choose A2 = A3 = 1, which is of course a possible 
choice. The rhs is bounded by 2, thus for the Ihs not to diverge we must also have (say) A4 = 1; but then the 
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State 




Sp 


5+ 


\W) = |011) + |101) + |110) 


1.516 


1.523 


1.527 


cosf |000) + sinf IVF) 


1.669 


1.669 


1.68 


cos^ ||000) +cosf |001) +sin||lll> 


1.425 


1.431 


1.431 


|0)(|00) + |11)) 


1 


V2 


V2 



Table 1: Sp and the bounds 5+ and Sp for some three-qubit states (states not- normalized). 

condition (j^) forces Ai = 1 too. In conclusion, if two eigenvalues are equal to 1, all the eigenvalues must 
be equal to 1. 

Since such constraints arc not easy to handle, it is interesting to provide a lower bound on Sp. A non-trivial 
one is obtained by simply choosing one possible realization of the eigenvalues, namely Ai = ±2, which due 
to ( p3| ) implies A2 = A3 = A4 = 0. This gives 

Sp ^ Sp — — max \txxx ^3:yy ^yxy ^yy2;| ■ (^^) 

^ {^,y} 

In most cases, we still have to rely on a computer program to calculate bounds ( p^ and (^5|). For these 
bounds, the optimization bears on nine parameters (for each qubit, two parameters define the (x, y) plane 
and a third one fixes the axes in the plane); while a direct optimization of (B3) bears on twelve parameters 
(two unit vectors per qubit). 

Let's conclude by a discussion of the quality of the bounds ( p^ ) and (p5|), based on some examples. Consider 
first the family of states |V'3(0)) (0)- ^ = 1 to « ^ we find S^ = S = S+ = 2 sin 20. For 

smaller values of 4>: (i) S^ = max(2 sin 2(/), cos 2(/)); that is, S^ corresponds to S^^^ defined in (pT|). (ii) As 
we discussed above, the exact value S essentially follows S^ , but for small deviations, (iii) The upper bound 
5"+ increases again, from S^{^) ~ 1.175 up to 5*^(0) = \/2. Thus this bound turns out to be too rough in 
the region (/> < where \ips{(j))) ceases to violate Mermin's inequality. 

In table 1^ we give Sp and the bounds 5+ and Sp for some other states. All the possible cases S^ = S < 5+, 
S~ < S = S^ and S~ < S < S^ are present. On all these examples, we see that at least one bound is very 
close, if not identical, to the exact value. 

4 Conclusion 

Bell's inequalities for systems of more than two qubits are the object of renewed interest, motivated by the 
fact that entanglement between more than two quantum systems is becoming experimentally feasible. A 
link between Bell's inequalities and the security of quantum communication protocols has also been stressed 
recently Q. 

Here we focused on inequalities obtained by measuring two observables per qubit, and we gave the spectral 
decomposition of the corresponding operators. With this tools, we studied the violation of Mermin-Klyshko 
inequalities for some states that exhibit N-qubit entanglement. We proved numerically for N = 3,4,5, and 
we conjectured for all A^, that there exist pure entangled states that do not violate these inequalities. 
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Appendix A: Relationships between Bn and B'^ 

We derive in this Appendix some properties of the MK operators Bn and B'^ that were not discussed in 
previous publications Lemma ^ was demonstrated independently, and with different mathematical 

tools, in H. 

Lemma 2 B^ — B'^^ for all n. 
Proof: From we have 

Bl - i(l + a„.a:,)l®i?2_i + ^(l-a„.a;)l®i3'Li - ^a,„Aa;® . (26) 

B'^ is obtained by exchanging the primed with the non-primed objects. Therefore 
Bl - B'J = a„ • a'^ (S^^ - S'Li) « {B^ ~ B',') = 
since it can be calculated explicitly that i?| = £?2^ = 1 + Ua^Aa'^ ® o'aiAa'^ ■ o 

Lemma 3 The explicit expressions for B^ and for the commutator [Bn , B'n] are given respectively by M. 



and ^(\)- The anticommutator is {i3„,i3^} = 2(a„ •a^)...(ai -a'j^) 1. As a corollary, note that Tr{BD — 2". 



Proof: Lemma H allows us to rewrite (p6|) as 



Bl = Bl^, ~ aa„Aa'„ [Bn-1 , SLl] • (27) 

Another standard calculation from (||) leads to 

[Bn , B,'J = 1 (8 [Bn-1 , + 2iaa^Aa',^ ® S^^, . (28) 

The structure of these two equations can be best seen by introducing the notations Bl = Pn, [Bn , B'n] = 
2iQn and cra„Aa^ = S„. We have then 

Pn = 1®P„_1 +S„®Q„_1 
Qn = l«)0„_i+S„®P„_i 

The recursive solution is a matter of patience. Using P2 = 1 + S2 Si and = 1 ® Si + E2 ^ 1 we find 

Bl = B'n = I2" + ^ CTa^Aa^CTa^Aa^. + ^ a^^a'.<^aj Aa'.^r Aa'i^CFaiAa[ + ■■■ (29) 
i<j i<j<k<l 

[Bn ,B'n] = 2t\J2 (30) 
\ i i<j<k J 

where the dots indicate the sums over all products of an even, resp. an odd, number of aa^Aa'.- 
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Finally, the anticommutator is also found through a direct calculation from (^: 



The conclusion follows from {i?2, -82} = 2(a2 • a2)(ai • a'j^) I4. o 



Appendix B: Some results about the maximal eigenvalue 

In the main text we exhibited a set of eigenvectors of B„ that have a remarkable symmetry. But this set 
would lose much of its interest if all eigenvalues were degenerate. Here we show that at least in one case 
(which is an interesting one) we can be sure that there are non-degenerate eigenvalues. 
Let's first sort the eigenvalues of i?„ in decreasing order: Ai > ... > By virtue of lemma = 

— In particular, Tr(i3^J = + A2 + ... + On the other side, we have noticed in lemma ^ 

that Tt{BI) = 2", whence 

A2 + a2 + ... + a2„_i =2"-^ (31) 

In particular, if Ai = 2~2-^ then all the other eigenvalues (except of course — Ai) are zero. In general, Ai > 1, 
where the equality holds only in the "classical" case Ai = ... = A2n-i = — A2ii-i+i = ... = — A2" = 1. Equality 
( ^ implies 

A^ + A^ < 2"-\ (32) 

whence the following: 

Lemma 4 let Ai and A2 he the two greatest eigenvalues o/_B„. If \i > 2^^^, then it is non- degenerate, and 
moreover A2 < 2^~^ . 

Now, if p is a n-qubit state exhibiting m-party entanglement, m < n, it can be shown that {Bn)p < 2^™^^)/^ 
[0 |l3|. Thus Ai > 2^^^ for _B„ means that we have a n-qubit violation. So our last lemma reads: if 
the parameters a of i?„ are such that a n-qubit violation is possible, the maximal eigenvalue of _B„ is 
non-degenerate. Actually we can prove even more: 

Lemma 5 If Xi > 2^^^, one cannot find two orthogonal states that both satisfy the condition 

(5„)^>2t-i. (33) 
Due to lemma^, the same holds for the condition {Bn)ip < ~22~^. 

To prove this lemma, we determine a necessary conditions for a state ji/i) to satisfy (p3[). Let's decompose 
\ip) on the basis of the eigenvectors of i3„: \tp) — VPi "^here is an eigenvector of i?„ for the 

eigenvalue A^. With these notations 

2" 

PiXi < piAi + (1 - pi)A2 . 
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Therefore if piXi + (1 — pi)A2 < 2"/^ ^, the requirement ( ^3|) cannot be satisfied. In other terms, a necessary 
condition for ( ^3| ) to be satisfied is 

Pi > ^ — = P (34) 

Al - A2 /il - fi2 

(we introduced the notation 2"/^^^/ii = Aj in order to show that p does not depend explicitly on the number 
of qubits n). It can be shown using ( p2[ ) that -i < p < 1. The limiting case p = 1 corresponds to Ai ~ 2"/^~^, 
in which case of course (^) cannot be satisfied. The roughest criterion that we can state is therefore the 
following: given _B„ such that Ai > 2"/^^^, a state \^) cannot satisfy ( |33| ) if pi = < i. This 

criterion is enough to conclude the proof of lemma ^. 
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